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Abstract 

We discuss the inclusive dilepton cross section for proton (quark)-nucleus colli- 
sions at high energies in the very forward rapidity region. Starting from the cal- 
culation in the quasi-classical approximation, we include low-x evolution effects in 
the nucleus and predict leading twist shadowing together with anomalous scaling 
behaviour. 



1 Introduction 



There is increasing evidence that hard probes PP are an excellent tool for analyzing the 
matter produced in high-energy heavy ion collisions at RHIC [2], especially when cali- 
brated against similar probes in proton-proton and proton (deuteron)-ion reactions [Hj. 
At central rapidities the fact that high-p^ hadron production in AA reactions is sup- 
pressed as compared to the production in pp collisions times the expected number of hard 
collisions, N co u, gave strong support to the idea that dense, hot matter is produced in AA 
collisions causing jets to loose a significant amount of their energy, while passing through 
the dense matter and before producing the observed high-p^ hadron. This picture was 
further confirmed when, at central rapidities, high-p^ hadron production in dA collisions 
did not show any suppression as compared to the expectation from pp. The lack of a 
suppression in dA as compared to pp collisions of course also means that there is little or 
no nuclear shadowing, at the hard scale determined by high-p^ hadron production, in the 
central rapidity region ^j. 

Recent dA data on high-p^ hadron production at large rapidity (toward the deuteron 
side) from the BRAHMS Collaboration [5] show a significant suppression of hadron pro- 
duction in dA collisions compared to the expectation from pp collisions. This result has 
aroused a lot of interest, because it suggests that there may be a significant amount of 
(leading twist) gluon shadowing in nuclear wavefunctions in the region probed by forward 
hard scattering at RHIC. The strong interest is connected to the fact that strong (leading 
twist) gluon shadowing appears difficult to understand outside of pictures which have 
gluon saturation H] (color glass condensate [HI M GDI )> which to a significant extent is 
driven by BFKL evolution [TT] . 

In many ways hard photon or /i-pairs coming from virtual photonsfHEIl are a better 
probe than high-p ± hadrons [HI [ISl [H3 [IZ1 [IHl ESI EDI EH E21 EI3 EH EHJ With hard 
photons (2HI EH EE1 EH [Ml EH one is less sensitive to fragmentation effects and final 
state effects are absent. This means that at transverse momenta around 2 — 3 GeV, 
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one can expect leading twist factorization to be accurate, and hence x-values of the 
gluon distribution of the nucleus down to values somewhat smaller than 10 -3 should be 
accessible. The main purpose of this paper is to explore, and estimate, the size of the 
suppression one might expect to see in such reactions. Our discussion is based on a picture, 
where the McLerran-Venugopalan model jHj is taken to represent the gluon distribution 
in a hard RHIC reaction at central values of rapidity, y = 0, and BFKL evolution [TT|l3*!2] 
is used to evolve the distribution to higher values of y. 

There has already been quite a lot of work studying hard photon and /i-pair production 
in p(d)A collisions j2H]- In a pioneering series of papers Kopeliovich and collaborators 
[23 12H| have studied Drell-Yan production in the RHIC and LHC kinematic regions using 
a dipole picture of the /i-pair production. Gelis and Jalilian-Marian [23 EOJ arrived at 
equivalent results in a color glass condensate picture, where the dipole of Kopeliovich 
et al. is replaced by a product of two Wilson lines evaluated in the field of the color 
glass condensate. Jalilian-Marian [25 then calculated the suppression factor, however, 
for the fc_i_-integrated yields in /z-pair production in p(d)A versus pp reactions taking the 
dipole cross section as determined by Iancu, Itakura and Munier [22] in fits to HERA 
data, based on the geometric scaling following from BFKL dynamics not too far from the 
saturation boundary of the color glass condensate. Quite a strong suppression is found in 
the analysis of [25 > because the gluon distribution used there has leading twist shadowing 
in contrast to the models in [23 EH 123 ED] 

In this paper we evaluate direct photon and /x-production in terms of standard fac- 
torization formulae. We remind the reader, how /c^-factorization formulae arise, and why 
fc_L-factorization is more efficient than ordinary operator product factorization, when one 
is dealing with small- x processes. Our general discussion is not tied to saturation or color 
glass condensate assumptions, but rather is a general leading twist discussion. However, 
because it is leading twist only, in contrast to previous discussions, it should only be used 
for moderate transverse momentum, say k\ > 4 GeV 2 . When we take the unintegrated 
gluon distribution, which appears in our formulation to be given in terms of the anoma- 
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lous scaling, which occurs in the BFKL based saturation picture jH2], our overall picture 
is very close to |3*T| . 

The outline of our paper is as follows: 

In Sec. 2, we derive a fcj_-factorized formula for high-fc^ transversely polarized (virtual) 
photons produced in a quark-nucleus (hadron) collision. The corresponding formula for 
lepton-pair production, with lepton pair mass M, is given by 

da gA^l+l-X q em do* A -~<* X 

d 2 bd 2 k ± d In zdM 2 ~ 3irM 2 d 2 bd 2 k L dhiz ) ^' 

where kj_ is the transverse momentum of the 7*, and z is the longitudinal momentum 
fraction of the 7* with respect to the incident quark momentum, z = k + /p + , where we 
have the limit p + — > 00 in mind, b denotes the impact parameter of the qA collision. 
Throughout the paper we shall refer to direct photon production, but lepton-pair produc- 
tion formualae follow easily from For simplicity we consider incident quarks rather 
than protons. 

In Sec. 3, we review the form that the unintegrated gluon distribution takes in sat- 
uration (color glass condensate) models. We do this first in the McLerran-Venugopalan 
model jH], which has gluon saturation but no gluon shadowing, and then for the case, 
where a significant amount of BFKL evolution is added to the McLerran-Venugopalan 
model, which is taken as the initial condition for that evolution. With BFKL evolution 
jTTJ |H2] gluon shadowing appears and the fixed impact parameter unintegrated gluon 
distribution scales with A (roue hly) like A~z~ , with A ~ 0.372. 

In Sec. 4, we present numerical results, which suggest a significant suppression of hard 
photons in the forward rapidity region in p(d)A collisions as compared to pp collisions. 
Our results, we hope, are encouraging for experimenters trying to measure the suppression 
at RHIC. 

In Appendix A we relate the k±- to the impact parameter representation. 
In Appendix B and C we relate our /cj_-factorized formulation to the more standard 
QCD factorization. We show explicitly that the anomalous scaling formulae, which ap- 
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pear in the fcj_-factorization formalism, lead to an (integrated) gluon distribution which 
obeys the renormalization group equation with an anomalous dimension given by BFKL 
evolution. 



2 Dilepton production cross section 

2.1 Factorized formula for the inclusive 7* cross section 

In the following we only consider the case for virtual photons with transverse polarizations, 
A = 1,2. The qA — > 7*X cross section is obtained from the /c^-factorized formula con- 
taining the photon phase space, the photon emission amplitude A\ and the unintegrated 
gluon distribution <f>a, explicitly 

where 2N~ ^^^ may be viewed as the differential high energy qA — > qA cross section. 
In the definition of (fro in © the integration with respect to the impact parameter is 
implied, 

</> G (q±,Y) = J d 2 b<f )G (b,q±,Y = hxl/x). (3) 

We shall discuss in the following subsection, how fcj_-factorization applies to qA — > 7*X 
at large k±, and how it leads to (j2j) (see also [HI]). 

The photon emission amplitude A\ is expressed in terms of the polarization vector e* ± A 
and the transverse momenta k±, and k{_ = k± — zq±, respectively as 

where we use the short hand notation [33] 

e ± A . {k ± } = (1 - z/2)e ± A • k± - iz/2 e ± A A k±, (5) 

and 

rf = (1 - z)M 2 . (6) 
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Performing the polarization sum we note that 

l + (l-z) 2 



EM 

A 



\ x ■ k\ 



k! 



k\ +r] 2 k'l + 



(7) 



Finally the transverse 7* production cross section becomes 
da qA^*x a m a„ r . .0, fd 2 q± 



d 2 bd 2 k 1 din. 



x 



z 2 q± 



[1 + (1 



(kl + V 2 )[(k ± -zq ± ) 2 + V 2 \ [ki + v 2 (k 



1 



a(b,q±,Y) 



1 2 



zgl) 2 + ?7 2 _ 



We remark that in the paper j2H] the longitudinal contribution is included 1 



d 2 bd 2 k 1 din z 



as -2(l-z) m±<f> G (bA±,Y) 



TlN, 

2 



XT] 



(9) 



_k{+V 2 {k ± -zq ± ) 2 + V 
In Appendix A we express the /c^-factorized cross section (jH} in terms of the dipole 
formulation in the impact parameter representation (2H1 EH EE] and we shortly mention 
the relation to DIS. Appendix B discusses the relation of (JBJ to the large kj_ LO pQCD 
cross section for the case of real (isolated) photons. 

2.2 Hard reactions and kj_— factorisation 

Before we continue and discuss the evolution of <pc with respect to increasing rapidity 
Y we briefly recount the origin and role of /c^-factorization in small-x hard reactions, 
especially the validity of Eq.©. 

In the usual QCD factorization jHEJ involving local gauge invariant operators in the 
operator product expansion at small- a; it may be necessary to resum a s In 1/x terms in 
both the coefficient functions and in the matrix elements, evaluated at a hard scale Q 2 , 
which occur. There is an alternative procedure in which the hard part of the reaction 
can be taken at lowest order in perturbation theory and the resummation done on what 
remains. In this k 1 -factorization a convolution in transverse momentum then remains 



1 The function C(l±) introduced in corresponds to °' — ~ J 
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to be done between the "factorized" parts while there is no convolution in longitudinal 
momentum because the formalism only exists in a leading order formulation. Indeed one 
of the shortcomings of the fc^-factorization formalism is that it is not known whether this 
leading order procedure is part of a more systematic procedure or not. On the other hand 
fc_i_-factorization is very useful when extremely small values of x are being considered where 
resummations in a s \al/x are paramount, which resummations are somewhat awkward 
in the standard hard QCD factorization plo] . 

Let us now examine how /c^-factorization comes about in the process of interest here, 
direct photon production in, say, quark-nucleon (or nucleus) collisions. The process is 
illustrated for a "typical" graph in Fig. 1. 
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— Y~ ~~? 7^ 
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Figure 1: Typical graph for photon production in quark-nucleon scattering. 

p is the incoming quark, P the target and k the hard photon setting the hard scale for 
the process. The lines q x , ...q r are gluons exchanged in the amplitude while q r +i, ■■■q n are 
the ones in the complex conjugate amplitude. We suppose that p + and k + are large, with 
z = k + /p + fixed, and we further suppose that all the gluons and quarks in the lower 
"blob" of the graph have + components of the momentum much less than k + . This latter 
assumption is important in /c_i_-factorization; a strong ordering in longitudinal momentum 
is necessary. In addition we suppose that k± is large while, for simplicity of discussion we 
take p± = 0. The lines, qi, which connect the hard part of the graph with the target, P, 
in general have q i+ k + ,p + . Now we limit our discussion to leading twist, in k\. In a 
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covariant gauge there may be many g^-gluons present, however in light cone gauge, with 
A_ = 0, the leading twist contribution can involve only two exchanged gluons [T31 
Thus, taking = we consider the graph shown in Fig. 2. 



p 




Figure 2: Two gluon exchange graph. 



There are three other graphs in the two gluon exchange or leading twist limit, however, 
they may be ignored as our object here is to explain /c^-factorization not to give a detailed 
calculation of terms in the factorized formula. 
In = light cone gauge the propagator is 



with rj ■ v = v_ for any vector v^. When applied to the hard part of the graph shown in 
Fig. 2, or to any of the other possible hard parts that may occur, it becomes [37| 



and similarly for Dppr(q), with rj ■ v = v + . The q a term in (jl(Jj) gives zero by current 
conservation, while fj a projects a relatively small component of the momenta. This leaves 
only the g + _ term in (JTUJ), which corresponds to (fTTf . 

While looks like a covariant gauge result this is not quite the case. In covariant 
gauge the leading twist contribution is not limited to the two gluon exchange term shown 
in Fig. 2; the many gluon exchange terms of Fig. 1 are also important. In addition, in the 



D aa '(q) — — g aa > — 



Vaqa' + Va'qa 



(10) 



v ■ q 




(ii) 



s 



present case, the lower blob in Fig. 2 must be evaluated in light cone gauge. A covariant 
gauge evaluation will give an incorrect result. 

Taking the graph shown in Fig. 2 along with the three graphs where either, or both, of 
the exchanged g-lines hook into the hard part of the graph before the photon, k, is emitted 
leads to the cross section formula given in Sec. 2. The unitegrated gluon distribution 
4>g{QLi y) is given by Va ' q V i' acting on the graphs in the lower blob of Fig. 2, as we 
now illustrate in Fig. 3. 




Figure 3: Unintegrated gluon distribution. 
Our normalization is such that 

(jfa'iq^Y) = — — -j (12) 

for a quark at lowest order in a s [2E]> and an additional factor N c for three quarks in a 
nucleon gives finally, e.g. (jB.6|) . 

For large Q 2 the following relation, using the qA — > qA cross section, between (fro and 
the gluon structure function may be derived 



Jo dQ±q ^ = 2N c J Mq± > Y)dq± 

(13) 

which, however, cannot be used in the scaling region. 

In contrast to normal hard QCD factorization fc^-factorization requires a convolution 
in transverse momentum be taken between the hard part and the unintegrated gluon 
distribution to arrive at a cross section, as given for example in (j2J). 
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3 Saturation and anomalous scaling 

In the leading twist region the function 4>a{b, q±,Y) in (jSJ), the unintegrated gluon struc- 
ture function 2 , is expressed in terms of the forward scattering amplitude N(b, x±,Y) = 
N q q(b, x±, Y) of a QCD qq dipole of transverse size x± with rapidity Y = In 1/x, scattering 
off a nucleus A at impact parameter b, by 



N r - - - 

^q ± ,Y) = —- ^- d 2 x ± e^V 2 x± N(b,x L ,Y) 



(14) 

N c 2 -* 2 f d 2 x± 



qlV 2 q J^e^N(b,x ±1 Y). 



(27T)3« S 

The function (j)c{b,q±,Y) may also be expressed in terms of the forward amplitude of a 
gluon-gluon dipole N G (b,x±,Y), obtained by replacing the number of colors N c by Cp 
gUllll- The relation fljfy is inverted by 

N &*^ = ^rj ^^fe.^-e^ 1 -^ 1 ], (15) 

illustrating the relation of 0g to the qA — > qA cross section more directly The phase 
factors in the bracket are due to the four graphs describing the different ways of gluon 
exchanges of the qq scatterings of the target nucleus |41j . 

Let us emphasize that the unintegrated gluon distribution (j>Q as defined in the previous 
section (see Fig. 3) is an object to be used only in leading twist fcx-factorized formulae. 
Eqs. (d3J) and (|T5j) are leading twist equations, valid in the scaling region - which we shall 
discuss later - and beyond. 

For illustration and later reference, we shortly summarize in Appendix B the pQCD 

leading order (LO) behaviour of N and 4>c- 

2 This function is sometimes denoted as modified gluon distribution h(q±,Y) |2()l . 
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3.1 McLerran-Venugopalan model 

Let us review the quasi-classical model by McLerran-Venugopalan jH] (at fixed b and at 
Y = 0) as a reasonable starting point of the Y evolution of 4>g, 

N MV (b,x ± ,Y = 0) = 1 -exp[-x 2 ± Q s 2 (b)/A}, (16) 

with the saturation scale |HZl EE] given by 

27r 2 (T 

Q s 2 (b) = ^ P T(b)(xG(x,l/xD), (17) 



with p the nuclear density and T(b) the profile function T(b) = 2y R\ — b 2 . xG is the 
gluon distribution in the nucleon. The unintegrated gluon distribution <f>Q V (b, q±,Y = 0), 
calculated from (fTIj). approaches (p G ° ()B.6|) at large q±, i.e. for q± 3> Q s , from above. 
The low momentum part is suppressed relative to the perturbative gluon; keeping Q s 
constant, independent on x±, one derives in this model for q± <^ Q s , 

w (g , f± , r = o) ^|. m 

As discussed in some detail in |2J , it may nevertheless be worthwhile to note that this 
behaviour differs from the one known from the unintegrated gluon distribution derived 
from the non-Abelian Weizsacker- Williams field of a nucleus, denoting it by <j) ww , which 
behaves in the quasi-classical approximation as 

cP WW (b, q u Y = 0) ~ J£- In (^i) 2 . (19) 

2n 2 a s q± 

Both functions, (fie and <p ww , do, however, agree for q± » Q s . When q± « Q s , all 
twists become important in this kinematic regime, and indeed, comparing (j!8)l and (j!9j) 
there is no unique prescription to define the unintegrated gluon distribution. In order 
to calculate the photon spectrum (J2J) at Y = we continue to use 4>g, as defined by 
(HU), together with ()16p. knowing that we will be interested in large enough transverse 
momentum only, so that in fact 4>g 
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The region around q± ~ Q s is enhanced, since the effect of multiple scatterings, 
resummed in f!16|) . rearrange the gluons in the nucleus ^mEHHO]- There is no shadowing 
in the quasi-classical approximation. 

Starting with (|16|) already at RHIC energies for dileptons produced in the central 
region (y 1 ~ 0) implies that the x- values in the gluon function in the nucleus are already 
small enough in order to justify the applicability of the McLerran-Venugopalan model, 
although there is no explicit x-dependence in this model. A rough estimate, based on 
hard two partons — > 7* /cj_-factorized kinematics, gives x ~ ^ at central rapidity; for 
collision energy y^s = 200 GeV and a transverse mass M± = 4 GeV a reasonable small 
value of x ~ 0.02 follows. This implies that the gluon number density at RHIC energies 
|2j is already large, i.e. saturated [Sid) 

/ 0^(6, gl, Y = 0)dV - (20) 

for fixed Q s (b). 

3.2 BFKL evolution in the presence of saturation 

Increasing the photon rapidity into the forward region, y 7 > 0, the values of x become 
rapidly small, namely x ~ {M±/ \/s)e~ yj , such that Y = lnl/x ~ y 7 increases with y 1 . 
In the following we fix Y = at y 7 = and treat Y as equivalent to y 7 for positive 
large rapidities. We work with the fixed coupling leading order approximation of the Y 
evolution, which effectively depends on the product of a s Y. 

In order to calculate the Y dependence of <pc we start from the BFKL evolution 
and write the amplitude N(b, x±, Y) in terms of the Mellin transform 

N(b, x ± ,Y) = -j — r(A - 1) exp [2a X (X)Y + (1 - A) \n(x 2 ± Q s 2 (b) / 4)] , (21) 

where we use the standard definitions at = a s N c /ir and the Lipatov function 

X(A) = ^(1) - ^(A) - - A), ^(A) = ^ . (22) 
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As usual the integration contour being parallel to the imaginary axis with < Re(X) < 1. 
Since we are in the following mainly interested in the region in which l/#±) is not 

very much larger than the saturation scale Q s , we keep in (|21j) the scale Q s (b) independent 
on x±. The normalization of (|21|) at Y = is given by the expression (|16|). This is best 
seen from the inverse Mellin transform in terms of the relation (t=x^_), 



oo 

A-2/i -t 







dtr-\i-e-*) = -r(A-i), (23) 



confirmed by partial integration. 

From the definition of (pa (|14J) and using 



dWe^Hxlr X = «^^(^) , (24) 

actually valid for 1/4 < Re(X) < 1, we obtain the Mellin representation of the uninte- 
grated gluon function, 

r(2-A)exo 2ay(X)Y — (1 — A) In I 

As a consistency check one obtains the result ()18j) at Y = by keeping only the pole at 
A = 2, dominating the behaviour at small values of q]_. Summing the contributions of all 
poles, A = 2, 3, one obtains 



A-l 



(25) 



N c ql \ q\ 



<p A J v (b, q ± , Y — 0) — 7 q \ exp 

27r a * Q 2 s (b) 



(26) 



Qs 2 (b). 

in case of a "frozen" scale Q s 2 . Inserting (|2T)j) into (JT^J) gives back (fTB^l . 

Up to the normalizing factors the function (ftc in (J2I3) has the structure of the amplitude 
T(Q,fi,Y) discussed in [H21I12], when identifying Q = q± and fi = Q s (p). 

Following the same steps as described in the paper [32], we consider the solution of 4>g 
for large values of a s Y extended into the geometric scaling region [33] • This is achieved 
by demanding that 4>G(b,q±,Y) vanishes close to the saturation boundary, i.e. for q\ < 
Q 2 s {b, Y), to be defined below. This pragmatic procedure includes non-linear effects which 
are present in the Balitsky-Kovchegov equation for N(b, x±,Y) [33]. This characteristic 
behaviour is also discussed in [33] from a more mathematical point of view. It is achieved 
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by a linear superposition of two BFKL type solutions by shifting the positions of their 
maxima by a finite amount. The final scaling solution, expressed in terms of the Y 
dependent saturation momentum 



Q 2 s (b,Y) = c s Q s \b) 



cxp 



2a ^Y 

1— Ao 



(a s Y) 



2(1-A ) 



(27) 



is (for the case of constant a s 
<f> G (b,q±,Y) ~- 



^max 



A ) exp 



111 



Ac) hi 



Q 2 s (b,Y)\ 



A, 



{21 



where c s and 0g ax = 0{l/a s ) are constants. The value of the anomalous dimension Aq is 



determined by 



x'(Ao) 



1 



-, A = 0.372. 



(29) 



x(A ) 1 - A 

Obviously 4>q is maximal, <pc = 0™ aX ' wnen Q± — Qs(b,Y), and 4>c = for q± < 
Q s (6,F)exp(-2( T rA^). 

It is well known that this leading order calculation with fixed coupling yields a large 
exponent in (|27j) . namely 2a = 4.66..a s , which is too large to agree with phenomenol- 
ogy jSSllini- However, this discrepancy is resolved in j!2], using the next-to- leading BFKL 
formalism, which as a result reduces the exponent to a value in agreement with the Golec- 
Biernat and Wiisthoff model [I'd] . 

It is important to note that this analytical function (f2*Kj) successfully compares with the 
numerical studies (201 HJj of the Kovchegov equation. Indeed in [20] a good fit by (|2*5jl is 
obtained for a fixed value of the anomalous dimension, Aq = 0.37, and for 5 < q±/Q s (Y) < 



1000, mainly because of the logarithmic factor, In 



i(b,Y) 



which is present in f!28j) . This 



comparison also indicates that the scaling behaviour is rather rapidly approached. 

The A, respectively the number of participants N part , dependence of the unintegrated 
gluon distribution (}2*Kj) is dominated by the behaviour for large A by 



<j> G (b = 0,q ± ,Y)~A-z SL \nA. 



(30) 
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rather than by A 1 / 3 . This is leading twist gluon shadowing due to the anomalous behaviour 
of 4>c with a non-vanishing value of Aq, in the extended geometrical scaling window 



Qt(b, Y)<q 2 ± < Qi(b, Y) exp ^a X "(X )Y , (31) 

and because of Q 2 s {b = 0, Y) ~ pRa — A 1 / 3 . We note that the upper bound given in (J3TJ, 
which has its origin in the diffusion region, is more stringent with respect to its Y depen- 
dence than the one quoted in g3], which behaves as Q*(b,Y)/Q 2 s (b) « exp [4a s ^Y]. 

Also when compared to the LO perturbative behaviour, <f) G ° ~ A given by (jB.lOj) . 
even stronger suppression of the gluon density (}3T)j) is observed. 

The consequences of these derived scaling properties of (f>a in (J2%|) for the dilepton 
cross section are analyzed in the next Section. 

4 Anomalous scaling and shadowing in dilepton pro- 
duction 

4.1 Qualitative results 

Before analyzing the transverse 7* differential cross section in the fc^-factorized form of 

(JBJ) in more detail we first investigate its scaling properties. We define 

da qA ^* x da qA ^* x 

<Pb =^ + r > 2) dVk ± dlnz • (32) 

Assuming r\ < k±, such that k± is the hard scale, we may approximate this cross section 

by 

d% = J ^ H (k±, *f±, z)<h$ ^/Qs(b, Y)) , (33) 

with 

H(k ± , q ± , z) = + (1 _ z f] S. , (34) 

(similar to the definition given in the previous section). 

Inserting the scaling function (J2EJ), the cross section (|33|) scales approximately as 
follows, 



d2b =Mk±/(zQ a (b,Y)),b)F 
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(35) 



where F is expected to be a slowly varying function of k±. In order to exhibit the 
anomalous A dependence we deduce a parametric estimate of the ratio with respect to 
the proton target, 

For central collisions, b = 0, and assuming that the extended geometric scaling regions 
for protons p and nuclei A indeed overlap, this ratio becomes 

R pA « A~ x °/ 3 . (37) 

Because of the nonvanishing anomalous dimension Ao, we thus predict shadowing of 7* 
production in quark-nucleus scattering at fixed k± and Y, at a constant level. The estimate 
(J37j) is based on approximating (J28|) by 

MkJQ.$Y)) « {kl/Q%Y)f°~\ (38) 

and on (c.f. (fT7|) and (j2T|) ) 

^(6,F)| A ^pT(6), (39) 

whereas the scale of the proton does not depend on b. 

A similar suppression in terms of anomalous scaling, as given e.g. by (|37[). is also 
predicted for the nuclear modification factor R^ A in case of gluon production [THl EE I2D1 

ED 1221 121 123- 

4.2 Quantitative results 

For illustration we present numerical estimates for the transverse 7* differential cross 
section (JHJ), for the RHIC energy y/s = 200 GeV. We are estimating the ratio R p a (}36|) as 
the ratio of central, 6 = 0, versus peripheral, b > 0, qA — > 7*X collisions, as follows 

da qA ^* x /d 2 b I da qA ^* x /d 2 b 



Rn 



M pT(b) centrai/ p T{b) 

where we choose, according to (JSHJ), 



(40) 

peripheral 



peripheral 



P T(b) Q 2 s (b = 0,Y) 

central 
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(41) 



The peripheral collision is assumed to be such, that N part = 1, i.e. the proton. For the 
numerics the ratio (J4*T]) is taken to be equal to {Q s ,pj Q^X^) — ^4 -1 ^ 3 - The central (periph- 
eral) 7* cross section in ()4L)p is evaluated at the scale Q s (b = 0,Y) (Q s (b,Y)\ periphera i = 
A~^Q s (b = 0,Y)). 




k T (GeV) 

Figure 4: McLerran-Venugopalan model: R p a as a function of k±. Solid curve for M = 2 
GeV and fixed z = 0.7; dashed (dot-dashed) curve forM = 4(2) GeV and y 7 = 3. 

In order to set the reference we give results based on the McLerran-Venugopalan model 
jH] as described in Sec. 3, when using (pQ V of (J2SJ). Instead of explicitly taking the scale 
(DID we fix the values at b = by QfJ = 1 GeV for A = 200, and Q S:P = 1/A 1 / 6 GeV 
~ 0.41 GeV for the proton target, respectively. 

For small photon rapidities, y 7 ~ 0, and dilepton masses of M = 2 — 4 GeV the ratio 
R p a of (j4T)j) as a function of kj_ is essentially R p a = 1. This is easy to see from (jHJ): small 
values of z < 0.05 imply that H of ()33|) is approximated by if ~ q\/ (k± + r] 2 ). Because of 
(JH it follows that Rp A w A-^^QfJ/Q,,,,) 2 ~ 1. For comparison with the prediction for 
evolved gluons we also consider larger values of y 7 , e.g. y 7 = 3, in this model. The results 
for R P A are plotted in Fig. 4 for two values of M = 2 and M = 4 GeV, respectively. A 
very small suppression, R p a > 0.95, is observed for < 1 — 2 GeV. Keeping, however, 
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1 . 



5 



. 1 1 1 1 1 1 1 1 1 1 

1 3 5 7 9 

k T (GeV) 

Figure 5: BFKL-saturation model: R p a as a function of k± for different values of y 1 and 
for M = 2 GeV: dotted (y 1 = 0.5), short-dashed (y y = 1.5), long-dashed {y 1 = 3.0). Solid 
curve for M = 4 GeV and y 7 = 3. 

the value of z fixed and large, e.g. z = 0.7, a significant Cronin type pattern for the 
transverse Drell-Yan spectrum is observed (solid curve in Fig. 4): shadowing for k± < 1 
GeV, and enhancement above. This confirms the results first presented in |27|. However, 
we have to keep in mind that this model is only reliable, when k± > Qg . 

In order to obtain results at large photon rapidities based on the BFKL evolution 
in the presence of saturation we pragmatically have to choose the Y dependence of the 
scale, instead of the one given by (|27|). As already discussed in Sec. 3 we take the one 
compatible with phenomenology, following |46j . 

Q 2 s (b = 0, Y = yj = (Q™ V ) 2 exp(\ GBwyi ), X GBW = 0.3 . (42) 

In accordance with the discussion in Sec. 5.1 significant shadowing is obtained, as 
shown in Fig. 5 as a function of k±, especially when the dileptons are produced rather 
forward, e.g. with y 7 = 3. Similar results, however, for fcj_-integrated dilepton rates are 
presented in j3*T| . 

When k± » Q s (b = 0, Y) the ratio R p a becomes essentially independent on the 
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5^ 



. 1 — 1 1 1 1 ' 1 ' ' ' 1 1 1 ' 1 1 1 

50 100 150 200 

A 

Figure 6: BFKL-saturation model: R p a as a function of A for = 5 GeV, y 1 = 3 and 
M = 2 GeV (solid curve). Dashed curve corresponds to A~ x °^ 3 . 

transverse momentum. 

Finally, we numerically check the statement in (|37j) concerning the A dependence of 
the 7* cross section ratio for central collsions. An example is plotted in Fig. 6, indicating 
that indeed the anomalous A dependence is to be expected for large nuclei, with the 
consequence of strong (leading twist) shadowing of photons/dileptons, when produced in 
the forward direction of pA, or dA collisions. 
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Appendix A 

In this appendix we relate the /^-representation of the 7* production cross section, Eq.(JHJ), 
to its impact parameter representation [SEIIZZIEH]- In order to go from k± to the conjugate 
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coordinate x±, we introduce the Fourier transform of the propagator 
e , A ■ k _l . f d 2 x± 



—i 



2tt 

'2 



K (r]x_ 



fci + if 

= —in I - 

2ir x± 

with x± = \x±\ and Ki the Bessel functions of the second kind, Ki(z) = — 4-Kq{z). 
The square of the radiation amplitude (JU) is expressed by 



(A.l) 
(A.2) 



E 



e , A • k i 



£; ± 2 + r/ 2 fc ± ' 2 + r/ 2 



(2tt) 2 
1-e 



K 1 (r]x ± )K 1 (r]y ± ) 



(A.3) 



+ 1 -e 



+*z<?X'?7x 



1-e 



-izgj.-(^x-yx) 



Following j27j, we introduce the transverse photon wave function for q — > 7*g, ip^ q {z, x±,r)), 
with 



A 



[1 + (1 _ z f] ^J^K^x^K^) . 



2-k* ■ x±y± 
Inserting (fTH|) for the unintegrated gluon function <pc, namely 



(A.4) 



d 2 q± 



(b,q ± ,Y) (1-e 



■N{b,zx ± ,Y) 



(A.5) 



q± 2na s 

into the 7* cross section (jBJ) we finally obtain, together with ()A.3j) the impact parameter 
representation of the 7* cross section [2*7] . 



D ik ± -(x ± -y ± ) 



(A.6) 



d 2 bd 2 k ± d\nz J (2tt) 2 
x ]T ^ 0, fx, £l, »?) \N(b, zS±, Y) + jV(6, zy ±1 Y) - N(b, z(x ± - y ± ), Y) 



In the gg-dipole approximation for scattering off a nucleus A, the amplitude N — 1 — S, 
where 

S(b,x ± ,Y) = exp[-^a q ,(b,x ± ,Y)} = exp[~x 2 Q 2 (b)/4} , (A.7) 
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is expressed by the a qq - cross section or by the saturation scale Q s 2 (b) = jf-Qs(b) given in 

(HZD [HZl EH1 

Let us add two useful relations: 



i) The limiting case for real photon emission, 77 — ^ 0, is obtained from (jA.4|) by 

£<&*±)<(*, ft.) - + (1 - zf\%^. (A.8) 

a Z7r x ±y± 

ii) Via crossing from q — > q r )*(M) into J*(Q) — > qq the DIS cross section is obtained 
from the expression (|A.6|) by the substitutions 

1 l + (l-Z) 2 x9 9 , A s 

, i.e. ^ ^ l-a 2 + a 2 , A.9 

1 — a z z 

and 

A'i(r/x ± ) = K x (y/(l - z)M 2 x 2 ^ -> #i (^/a(l - a)Q 2 xi) , (A. 10) 

where the gg separation is given by x± = = zx±, such that 

N(b,zx±,Y) -> JV(6,xj.,y) . (A.ll) 

The explicit expressions for deep inelastic scattering may be found in [2o1 137j . 

Appendix B 

Here we give some details on the behaviour of the leading order (LO) cross section for real 
and isolated photons produced via the process qA — > 'jX j^DJ. The photon takes away a 
large transverse momentum k±, k± 3> Q s , opposite to a recoil quark jet. 

The differential cross section for real photons, M 2 = 0, in the /c^-factorized form reads 
(see Eq. ©), 

da _ a em a s 1 + (1 - zf r 2 (j) G (b,q±,Y) , R ^ 

d 2 bd 2 k ± d\nz nN c k\ J 9± [q ± -k ± /z} 2 ' ' 

Assuming k± 3> zq±, i.e. such that collinear quark-isolated photon configurations are 

suppressed, we may write (|B.1|) as 



da OL em a s z 2 [l + (1 — z] 
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.r 



(PA^dlnz k\ \ z(l — z 
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xGJi= Q 2 = fci , (B.2) 



where xGa is the gluon distribution in the nucleus A, 

xG A (x, kl) = J d 2 b<p G (b, q±, Y). (B.3) 

in agreement with (|13|). For real photons xt = k±/y/s and z = xj_exp?/ 7 , such that the 
x value in the gluon is decreasing with increasing photon rapidity y 1 . The expression 
(jB.2|) is the same as derived from the hard LO pQCD 2 — > 2 Compton process qG — > 7g 
for production of isolated photons at large k± [HE]: indeed the corresponding hard cross 
section is given by 



da qG ~* iq na err ,ote 1 



em^s 

2 



S t 

i s 



dt N c s 

na em a s z 2 (l - z)[l + (1 - z) 2 } 

= ~N C ¥ x ' (R4) 

since s = z ^ z y i = —k\jz. After folding (jB.4|) into the expression for qA — > 7X with 
the help of the gluon function xG a in the nucleus A and performing the integration over 
the recoiling quark leads to (jB.2D . 

Using <mj), we relate the gluon distribution (frc to the dipole amplitude N q q in the 
quasi-classical approximation [S]. With (|A.7j) one obtains, 

N r - 

<p G (b,q±,Y) = — -f- J d 2 x ± e^V 2 x± [l-exp(-x 2 ± Q 2 s /4)}, (B.5) 
with the scale given by (fT7j) . At very large q± ^> Q s one finds, 

^ G °(b,q ± ,Y) = pT(b)^^. (B.6) 
This is derived by using the LO gluon distribution, 

xG LO (x, 1/4) = In * (B.7) 

7T X 1 ./V 



and thus 

N LO (b,x ± ,Y) =x 2 ± Q 2 (b)/A = %i p T(6)xiln-i^. (B.8) 



]6j) follows by using |1S] 



d x± jff 



ln ZTUT = —T' (B.9) 



(2vr) 2 V^A 2 y yrgi 
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valid for gl 7^ 0. Since 

/ &°(b, q ± , Y)d% = A^^\, (B.10) 
J 71 q ± 

the LO factorized result for the gluon in nucleus A, 

xG™{x, kl) ~ A {N \ 1)as In % = AxG LO (x, k 2 ± ), (B.ll) 

Z7T A^ 

in terms of the LO gluon in the nucleon, is finally obtained. 

Appendix C 

In this appendix we extensively discuss in more general terms the relationship between 
the fc_|_-factorization, which has been used in the previous sections and the usual QCD 
factorization [3B] involving local gauge invariant operators which appear in a Wilson 
operator product expansion. As in Appendix B, where we investigate this relationship 
for large k± » Q s in LO pQCD, we concentrate explicitly on the case of real photon 
production. 

C.l The forms of factorization 

We restrict here our discussion to the scaling region where the hard scale is above, but not 
too far above, the saturation region. When the hard scale is below the saturation momen- 
tum neither fc^-factorization nor the usual QCD factorization is applicable as higher twist 
terms are coherent with leading twist terms and all terms must be considered together. 
While it does make sense to talk of a gluon distribution which has reached saturation, 
that distribution does not appear simply in factorization formulae. When the hard scale 
is very large and outside the scaling region fc^-factorization may still be useful, but the 
issues are more straightforward than in the scaling region. 

We write generically a dimensionless observable in the fc^-factorized form 

a(Q, Y) = [ H(q ± , Q) <f> G (Q±, Y) , (C.l) 
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and in the QCD factorization form 

a(Q, Y) = ^j™ dy H(y) xG(Q, Y - y) , (C.2) 

where Q is the hard scale of the reaction, H is the hard part in the /cj_-factorized form, 
and H the hard part in the usual factorization. Here, and in the following, we suppress 
writting explicitly the dependence on the impact parameter b±. We suppose that the 
hard part in the fc_i_-factorized expression is not too nonlocal in rapidity while we cannot 
suppose such is the case for H. Finally the coupling in H and H should be taken at the 
hard scale Q. In our example of direct photon production the hard scale Q becomes the 
transverse momentum, k±, of the photon, while 

, {k± ,Y)^ {b± , k± ,Y) = d2 J (C.3) 



and (c.f. (TRT 



^[i + (i-*) V T / ^ (a4) 

iV c {q± - k±/zy 



In the scaling region we approximate 4>g and xG, respectively, from (|28j) by neglecting in 
the following possible constants under the logarithms. We write it in the form 

y ) = - (^P) ° H&/&.00) , (c.5) 

and 

xG(Q, Y) = ^ 1 A ° HQ 2 /Q 2 S (Y)) . (C.6) 

In this region the normalizing factors C and C are not necessarily the same, and actually 
we have not been able to relate them. 

C.2 The renormalization group 



In this section we show that (jC6|) obeys the renormalization group equation. In showing 
this we shall find a relation which will be crucial in relating H and H, which appear in 
(jC.ljl and (|C.2jl . Now in BFKL dynamics there are two alternate forms for xG, 

xG(Q,Y) = J ^ Ane ^nQ 2 /» 2 +(n-i)Y ; (C 7) 
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and (c.f. (jZEJ)) 

X G(Q, Y) — [ ^B x e^xWY+xmQVS ? (Q g) 

where the scale n is introduced to create dimensionless quantities, but xG does not depend 
on fi. The integral in (jC.7|) goes parallel to the imaginary axis with Re(n) to the right 
of all singularities of 7 n and A n in n. As discussed in Sec. 3, Eq. flC.8|) is, of course, 
not a perfectly correct representation of BFKL dynamics in the presence of saturation, 
nevertheless the aspects of f)C.8|) . which we use remain true even when the BFKL equation 
[TT] is replaced by the Kovchegov equation iTTl. 

In the scaling region the integrals in (|C.7|) and (|C.8J) are dominated by saddle points 
at n = n and A = A , where A satisfies (J2HJ). Also for (jC.7|) and (jC.8|) to describe the 
same function it must be true that 

n-l = 2ax(ln) , 
which determines 7„, while at the saddle points 

n - 1 = 2ax{\ ) , 

and 

Ao = 1n ■ 

The renormalization group equation jHEj is 

Q*—xG(Q,Y) = j o d yi {y)xG{Q,Y -y), 
with 7(1/) the gluon anomalous dimension 

Using (jC.6J) along with the result ()27|). 

[a s y]2(i-A 0) 
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(C.9) 

(CIO) 
(C.ll) 

(C.12) 
(C.13) 
(C.14) 



it is straight forward to get 

xG(Q, Y-y) = xG{Q, Y) e~ mxiXo)y 



1 + 



2a X (A ) 

i-Ao y 



(C.15) 



so long as y/Y <C 1. Using ()C.6J) on the left hand side of (|C.12jl . and ()C.15Jl on the right 
hand side one easily sees that (jC.12|) is satisfied if 



A = / dy e 
Jo 



-2ax(\o)y 



i{y) , 



and 



2ax(A ) 



dyyi{y) e~ mx{Xo)y 



(C.16) 



(C.17) 



1 - A Jo 

are both true. Ea. (IC~T6j) follows from (fHTOji . (fTTTT^ and (IC~T3l) . Ea. flHTH) can be written 

as 

i = - , (c.is) 

1 — An an 

which requires the inverse representation of (jC.13|) . The validity of (jC.18|) then follows 
from differentiating ()C.9J) with respect to n and using (I29|) . 

C.3 The relationship between /c^-factorization and QCD factor- 
ization 

We now reach the main topic of this section, namely the relationship between the two 
forms of factorization exhibited in (|C1|) and (jC.2|) . We begin with (jC.l|) and insert (jC.5|) 
for (j) G . Thus 



a(Q,Y) = 
Using ()C.6|) we arrive at 



— r H (<1±> Q) — — — 



Hil/Qt(Y)) 



(C.19) 



HqI/Q 2 ) 



MQ 2 /Q 2 S ( Y )) 

Now we may rewrite, using (|C.15J) , QCD factorization as given in (|C.2J) as 

2«x(Ao) 



(C.20) 



_, 1 roc _, 

<r(Q, Y) = xG(Q, Y)—j Q dyH(y) e -**(*0» 



l-Ao 



in(gvg2 (F)) 



(C.21) 
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Using (jC.lOJ) one can write (|C21J) as 



1 

Q 1 



with, of course, 



<t{Q,Y)=xG(Q,Y)- 



H n — 



re— 1 d 
1— Aq dn 



dyH{y) e~ (n ~ 1)y . 



n=no 



Comparing (fTT20jl and (fTT22jl it is easy to see that they are equivalent if 

C f d 2 q ± (Q 2 



H 



"0 



C 



irql \q 2 ± 



and 



n 



1 9 H 

no ~ C 



C rd 2 q± ( Q 2 ^ l X{) 



A dn 



{ \\ Hqi/Q 2 )H(q ±J Q) 



are satisfied. Eq. (jC.24|) is easy to satisfy, if we choose to define H(y) by 

dn 



H(y) 



2ni 



e {n - 1)y Hr, 



and with (fTTTTTf) 



(C.22) 



(C.23) 



(C.24) 



(C.25) 



(C.26) 



(C.27) 



Then it is straight forward to see that (jC.24|) is satisfied while ()C.25|) follows by using 



d"?n 



dn 



A ) 



(C.28) 



n=no 



n - 1 

which, after using (jC.10|) . is the same as (|0.18jl . 

Thus we see that fcj_-factorization as expressed in (|C.1J) leads to QCD factorization, as 
expressed in flC~2j) . with H(y) defined by (ITl26|) and (f(T27|l . While H(q ± ,Q) is a lowest 
order expression, H(y) is determined in terms of a resummation dictated by (|C.26|) and 
(IC.27J) and cannot be limited to its lowest order term. The simplicity of fcx-factorization 
is that all resummations are put into (fro with H remaining a relatively simple quantity, 
that is the hard part defining the reaction is more visible in /c^-factorization than in 
QCD factorization. We remark that including a common additional constant under the 
logarithms in (|C.5|) and (|C.6|) . respectively, does not destroy the derivation given above. 
It remains a challenge to understand how to extend /cj_-factorization beyond a leading 
logarithmic formalism. 
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